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Supernova pointing with low- and high-energy neutrino detectors 
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A future galactic SN can be located several hours before the optical explosion through the MeV- 
neutrino burst, exploiting the directionality of ^-e-scattering in a water Cherenkov detector such 
as Super-Kamiokande. We study the statistical efficiency of different methods for extracting the 
SN direction and identify a simple approach that is nearly optimal, yet independent of the exact 
SN neutrino spectra. We use this method to quantify the increase in the pointing accuracy by the 
addition of gadolinium to water, which tags neutrons from the inverse beta decay background. We 
also study the dependence of the pointing accuracy on neutrino mixing scenarios and initial spectra. 
We find that in the "worst case" scenario the pointing accuracy is 8° at 95% C.L. in the absence of 
tagging, which improves to 3° with a tagging efficiency of 95%. At a megaton detector, this accuracy 
can be as good as 0.6°. A TeV- neutrino burst is also expected to be emitted contemporaneously 
with the SN optical explosion, which may locate the SN to within a few tenths of a degree at a 
future km 2 high-energy neutrino telescope. If the SN is not seen in the electromagnetic spectrum, 
locating it in the sky through neutrinos is crucial for identifying the Earth matter effects on SN 
neutrino oscillations. 



I. INTRODUCTION 

Observing a galactic supernova (SN) is the holy grail of 
low-energy neutrino astronomy. The question "how well 
can one locate the SN in the sky by the neutrinos alone?" 
is important for two reasons. Firstly, the MeV-ncutrino 
burst precedes the optical explosion by several hours so 
that an early warning can be issued to the astronomical 
community Q, , specifying the direction to look for the 
explosion. Secondly, in the absence of any SN observation 
in the electromagnetic spectrum, a reasonably accurate 
location in the sky is crucial for determining the neutrino 
Earth-crossing path to various detectors since the Earth 
matter effects on SN neutrino oscillations may well hold 
the key to identifying the neutrino mass hierarchy 0, 0, 

liSl. 

Nearly contemporaneously with the optical explosion 
an outburst of TeV neutrinos is expected due to pion pro- 
duction by protons accelerated in the SN shock . This 
neutrino burst could produce of order 100 events in a fu- 
ture km 2 high-energy neutrino telescope, allowing for a 
pointing accuracy of a few tenths of a degree. Apart from 
the precise SN pointing, the detection of high-energy neu- 
trinos would be important as the first proof that SN rem- 
nants accelerate protons. 

The optical signal from a SN, if observed, can give the 
most accurate determination of its position in the sky. 
Apart from the observations at the optical telescopes, the 
multi-GeV to TeV photons associated with the acceler- 
ated protons in the shock could be detected on the ground 
by air Cherenkov telescopes after the SN environment be- 
comes transparent to high-energy photons. If a suitable 
x- or 7-ray satellite is in operation at that time, the SN 
would be visible in these wavebands starting from the op- 
tical explosion. A satellite like INTEGRAL could resolve 
the SN with an angular resolution of 12 arc- minutes [Icj . 

However, it is possible that the SN is not seen in the 
entire electromagnetic spectrum. This can be the case if 
it is optically obscured, no suitable x- or 7-ray satellite 



operates, and the air Cherenkov telescopes are blinded 
by daylight or the satellites and telescopes simply do not 
look in the right direction at the right time. It is also 
possible that not every stellar collapse produces an ex- 
plosion so that only neutrinos and perhaps gravity waves 
can be observed. In such a scenario, the best way to lo- 
cate a SN by its core-collapse neutrinos is through the 
directionality of ve~ — > ve~ elastic scattering in a water 
Cherenkov detector such as Super-Kamiokande 0, . 
Much less sensitive methods include the time-of-arrival 
triangulation with several detectors [ill IT^ | or the sys- 
tematic dislocation of neutrons in scintillation detectors 
that measure v e p — > ne + and the subsequent neutron 
capture [u| . 

The pointing accuracy of Super-Kamiokande or a fu- 
ture megaton detector such as Hyper-Kamiokande or 
UNO is strongly degraded by the inverse beta reactions 
D e p — ► ne + that are nearly isotropic and about 30- 
40 times more frequent than the directional scattering 
events. Recently it was proposed to add to the water 
a small amount of gadolinium, an efficient neutron ab- 
sorber, that would allow one to detect the neutrons and 
thus to tag the inverse beta reactions . Evidently this 
would greatly improve the pointing: a tagging efficiency 
of 90% would double the pointing accuracy UJ. At high 
tagging efficiency, however, the nearly isotropic oxygen 
reaction u e + 16 O — > X + e~ remains as the dominant 
background limiting the pointing accuracy. In this pa- 
per we analyze the realistic pointing accuracy of a water 
Cherenkov detector as a function of the neutron tagging 
efficiency. 

The directionality of the elastic scattering reaction is 
primarily limited by the angular resolution of the de- 
tector and to a lesser degree by the kinematical devi- 
ation of the final-state positron direction from the ini- 
tial neutrino. Extracting information from "directional 
data" is a field in its own right 0, 0|. An efficient 
method is the "brute force" maximum likelihood esti- 
mate of the electron events, taking into account the an- 
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gular resolution function of the detector on top of a nearly 
isotropic background. For a large number of events, the 
accuracy with this method in fact asymptotically ap- 
proaches the minimum variance as given by the Rao- 
Cramer bound (l8L fisj . However, for a small number of 
signal events, N s < 200, we find that a fraction of the in- 
formation content of the data as measured by the Fisher 
information |20j| cannot be extracted by even the maxi- 
mum likelihood method. Using the Rao-Cramer bound 
therefore overestimates the pointing accuracy of an ex- 
periment for small N s . In this paper, we determine the 
realistic accuracy by using a concrete and nearly optimal 
estimation method. 



Since the angular resolution depends on the event en- 
ergy, the likelihood method requires as input the func- 
tional form of the neutrino energy spectra that are only 
poorly known. It is difficult to systematically take into 
account the errors introduced by a wrong choice of the 
fit function for the likelihood method. Therefore we look 
for "parameter-free" methods like harmonic analysis that 
use only the information contained in the data and ex- 
ploit the symmetry of the physical situation. We dis- 
cuss the efficiency of two methods closely related to the 
harmonic analysis and find a simple iterative procedure 
making them nearly as efficient as the maximum likeli- 
hood approach. We use the most efficient method thus 
obtained for analyzing the simulated events at a detector. 
This makes our estimations realistic and even a bit con- 
servative, since the existence of a better parameter-free 
method is not excluded. We also study the dependence 
of the pointing accuracy on the neutrino mixing parame- 
ters and the initial neutrino spectra, and use the "worst 
case" scenario in order to estimate the accuracy. 



Finally, we briefly study the pointing accuracy in high- 
energy neutrino telescopes that can pick up the TeV neu- 
trino burst expected around the time of the optical explo- 
sion. While an accurate pointing is easy for any of the 
existing and future neutrino telescopes if the neutrinos 
are observed through the Earth, it is far more difficult 
against the atmospheric muon background from above. 
Even this would be possible for future km 2 detectors such 
as IceCube or Nemo that could detect around 100 SN 
events with TeV energies within about one hour. 



We begin in Sec. [H] with a discussion of the statistical 
methodology for extracting information from directional 
data using a toy model. In Sec. lIIII we study the realistic 
SN pointing accuracy of water Cherenkov detectors as a 
function of the neutron tagging efficiency, using realistic 
SN neutrino spectra. In Sec. II VI we turn to the pointing 
accuracy of high-energy neutrino telescopes. Sec. is 
given over to conclusions. 



II. ANALYZING DIRECTIONAL DATA 
A. Pointing with maximum likelihood estimate 

As a first case we study SN pointing with the v-e elastic 
scattering events in the absence of any other background, 
thus obtaining a bound on the realistic pointing accuracy. 
To this end we use the toy model introduced in Ref . , 
i.e. we imagine a directional signal that is distributed as 
a two-dimensional Gaussian on a sphere. This choice is 
motivated by the observation that the scatter of signal 
event directions is dominated by the angular resolution 
of the detector, and by the assumption that the angu- 
lar resolution function is Gaussian. Later in Sec. IIIII we 
consider a more realistic approximation to the angular 
detector response. 

The angular width of the assumed Gaussian distribu- 
tion is denoted by S s , where s stands for "signal." As a 
further simplification we assume 8 S <C tt/2, allowing us 
to approximate the sphere by a plane. Taking the sig- 
nal to be centered at i?o = 0, the probability distribution 
function (pdf) of the signal events is 

f a (0, 4>)ddd<t> = \ exp dfi (1) 



C 



2S'i 



with d^i = sinMM(f>. Here C = J dfj, cxp[-tf 2 /(2(5 2 )] is 
a normalization constant taking the value C = 27r<5 2 for 
planar geometry. 

In the case of Super-Kamiokande, around 300 elastic 
scattering events constitute the directional signal. As- 
suming the mean electron energy to be 11 MeV, a cone 
with opening angle ^68 ~ 25° around the true direction 
contains 68% of the reconstructed directions [2l|. Solving 



(£68*7180°) 



ddf s (d,4>) = 0.68 



(2) 



we find S s ~ 17°. For N s ^> 1 signal events with- 
out background, the Central Limit Theorem implies that 
the mean reconstructed direction is within S s /\/N s of 
the true direction for 68% of all SN obervations. This 
quantity, which is 1° for 300 events, gives the absolute 
lower bound on the pointing accuracy in the absence of 
all backgrounds. We note that Ref. uses $s ~ 25° 
and hence obtains 1.5° for the pointing accuracy. We 
further note that our 8 s /\/Nl ~ 1° implies that 95% of 
all SN reconstructions lie within a circle of angular radius 
2.4° of the true direction. 

The main degradation of the pointing accuracy is 
caused by the nearly isotropic inverse beta decay back- 
ground. The extent of this degradation was addressed 
for the first time in Ref. [llj . The pdf on a sphere that 
represents N s signal events distributed like a Gaussian 
around the direction ($Oj0o) a s well as the Nf, isotropic 
background events is 



dfi 



Nb + N s 



N b , N s 
4^ + ^ 6XP 
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(3) 
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where 

I = cos -1 [cost? cos $0 + sin $ sin $0 cos(0 — <f>o)] (4) 

is the angular distance between the direction of an incom- 
ing neutrino and the experimentally measured direction 
of the Chercnkov cone. We have introduced here the 
usual notation f(x\xo) for the pdf to stress the depen- 
dence of / on the data x = (</>, -d) and the parameters 

Xq = ((po^o). 

The maximum likelihood estimate (MLE) method is 
the most efficient way to extract information from statis- 
tical data. For the pdf of Eq. © the likelihood function 
for N events is 



N 



(5) 



where ($<- a \<j)^) are the coordinates of the a th event. 

One commonly uses the Fisher information matrix poj 
to estimate a lower bound on the uncertainty of the pa- 
rameters extracted by the MLE method. In our case the 
Fisher matrix is defined as 



F 



& 2 lnL^o^o) 
dGidQi 



(0) 



where i,j = 1,2, Q 1 = i} , 6 2 = (f>a and (...) denotes 
an average with respect to fd"dd<j). Since the off-diagonal 
elements of this matrix vanish, the error in the measure- 
ment of the two angles is simply given by 



A6, = ^T/E n 



(7) 



This lower bound on the pointing error is also known as 
the Rao-Cramer bound 0, 0] . 

For the sake of definiteness, we have chosen the SN 
direction to lie in the equatorial plane so that Ai?o = 
A(f>o, and define the pointing error At? as 



At? 



1 



N - 



(8) 



where $ is the estimate for z?o from a given method and 
N is the number of simulations used. The Rao-Cramer 
bound corresponds to the inequality 



(Atf) 2 > (A<?)| isher 



= 1/F, 



(9) 



where F may be F%i or F 2 2- 

The accuracy of the MLE method approaches asymp- 
totically the Rao-Cramer bound for a sufficiently large 
number of events |22l | . For a finite number of data points, 
the MLE saturates this bound only if the pdf used in the 
MLE can be written in the product form 



f(x\x ) = g(x )h(x)e-xp[A(x )B(x)} 



(10) 



If the pdf cannot be written in this form, the efficiency of 
the MLE is significantly smaller than unity at low values 
of N„. 



In Fig.QJwe show the efficiency e of the MLE in "Fisher 
units" , 



e = (A^ ishcr /(Atf) 2 



(11) 



as a function of the number N s of signal events while 
keeping the background-to-signal ratio /N s fixed at 30, 
which is the expected ratio of the inverse beta de- 
cay events to the elastic scattering events in a water 
Cherenkov detector. Though the MLE efficiency tends 
asymptotically to the Rao-Cramer bound for large val- 
ues of N s , this bound overestimates the MLE pointing 
accuracy by - 10% for N s < 200. 

We stress that the Rao-Cramer bound as an estimate 
of the MLE pointing efficiency is useful only in the limit 
where the spherical nature of our problem can be ap- 
proximated by a planar geometry, i.e. when the angular 
resolution of the detector is much better than 90° . This 
condition is satisfied in our case. For the general case 
of a pdf defined on a sphere the MLE still provides an 
optimal method to determine the pointing accuracy, but 
the Fisher information matrix no longer provides a direct 
asymptotic bound on the pointing accuracy. We are not 
aware of a generalized version of the Rao- Cramer bound 
that would relate directly to the variance of the pointing 
estimate in the case of a truly spherical problem. 




FIG. 1: Efficiencies of different estimation methods described 
in the text for Nt/N 3 — 30. RC corresponds to the Rao- 
Cramer minimum variance. 



B. Efficiencies of parameter- free methods 

The MLE is an optimal method to extract informa- 
tion from experimental data if the probability distribu- 



4 



tion function is known. This is not the case in our situ- 
ation, where the exact forms of the neutrino spectra are 
needed and these are only poorly known. It is therefore 
worthwhile to look for other methods which may be less 
efficient, but which do not depend on the exact form of 
the pdf. In our case of SN pointing we wish to consider 
methods that do not depend on prior knowledge of the 
exact neutrino energy spectra. 

Let us consider two pointing methods that exploit the 
symmetries of our physical situation, but are indepen- 
dent of the exact details of the pdf. If the efficiency of 
such a method turns out to be comparable to the maxi- 
mum likelihood method for the toy model of the previous 
section, then that method may be expected to be efficient 
also in the realistic case where the exact pdf is not known 
and the MLE method cannot be employed. 

One obvious approach is the "center of mass" (CM) 
method. The center of mass S of the events, 



N 

E 



(a) («) 



1,2,3 



(12) 



is taken to be the estimator of the true center of the dis- 
tribution, where x = (sin $ cos (f>, sin i) sin <f>, cos $) . For 
an ideal detector the event weights mS a > are equal and 
can be set to one. More realistically, the detection prob- 
ability p is a function of the detection angles <p and the 
weight is = l/p(tf( Q \ <^ Q )). 

A second approach is the "orientation matrix" (OM) 
method |l6j . Here, the major principal axis of the orien- 
tation matrix 



N 



M«) T («) 



(13) 



is taken to be the estimator. 

These two methods are equivalent to a harmonic anal- 
ysis, 



A? 



aw, 



(14) 



a=l 



restricted to the first moment for CM, and up to the 
second moment for OM. This equivalence can be seen 
either by a direct evaluation of Eq. Ijl4(l or by identifying 
Si with a dipole moment and relating Ty to a quadrupolc 
moment Qij . Since 3T^ = Qy + NSij , the direction of 
the major principal axis is identical for both ellipsoids. 
Note that the orientation matrix is a reducible tensor 
and therefore contains information from the first as well 
as the second moment. 

Neither of these methods requires any prior knowledge 
of the neutrino spectra or cross sections. However, they 
involve some loss of information and hence will give larger 
pointing errors than the MLE. In order to quantify the 
efficiency of these methods we generate a data sample 
according to the pdf of Eq. (|3J) and show the respective 
pointing errors in Fig. [21 Wc keep the number of signal 



events fixed at N s — 300, and show the pointing error 
A$ as a function of Nf,/N s . Note that in the absence 
of neutron tagging this ratio is expected to be around 
30-40. 




20 30 
Nb/Ns 



FIG. 2: Pointing error Ai9 for different estimation methods 
for N 3 — 300. RC corresponds to the minimum variance as 
given by the Rao-Cramer bound. 

Figure shows that the error of MLE is almost the 
same as the Rao-Cramer (RC) bound. However, the er- 
rors of CM and OM are much larger. One may also notice 
that OM is more accurate than CM. This difference may 
be attributed to the fact that whereas CM tries to exploit 
the spherical symmetry of the background, OM exploits 
the cylindrical symmetry of the background about the 
arrival direction, which is broken much more weakly by 
statistical fluctuations in the background. Moreover, in 
terms of a harmonic analysis, OM involves information 
from both I = 1 and 2 while CM involves only 1 = 1. 

In order to increase the efficiency of CM and OM, we 
use the physical input that the signal is concentrated 
within a small region around the peak. Cutting off the 
events beyond a certain angular radius would then in- 
crease the signal to background ratio and the above meth- 
ods may be applied iteratively to this new data. This 
procedure converges quickly and gives a much better es- 
timate of the incoming neutrino direction. The optimal 
value of the angular cut has a very weak dependence on 
the number of events and the background-to-signal ratio. 
It depends mainly on the value of i5 s and is found to lie 
between 2S S and 35 s . Within this range, the efficiency 
depends only weakly on the exact value of the angular 
cut. We tried both a sharp cutoff and a Gaussian weight 
function; both choices give practically identical results. 

The optimal value of the cut also increases slowly with 
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decreasing background-to-signal ratio, and in the limit 
of zero background, the method without cut is clearly 
more efficient than the method with cut since the latter 
now cuts off signal but no background. However, the 
variation due to changing the cut is of the order of only 
a few percent. Therefore, we keep the value of the cut 
to be constant at 40° for the analysis in this section. 
Our results for the pointing accuracy will therefore be 
somewhat conservative. 

We denote the CM and OM methods with this cut- 
ting procedure by CMc and OMc, respectively. Figure [21 
shows how the pointing error decreases drastically with 
the cutting procedure. It may also be observed that the 
accuracy of CMc stays close to that of MLE for low val- 
ues of the background-to-signal ratio, while the accuracy 
of OMc is close to MLE in the entire Nb/N s range. 

The efficiencies of all methods depend on both N s and 
Nb/N s . In Fig. Elwc show the efficiency e for different 
methods as a function of N s , keeping Nh/N s fixed at 30. 
For a large number of signal events, N s > 300, all meth- 
ods tend to their asymptotic efficiencies £oo. The OMc 
method is close to its asymptotic efficiency 0.77 al- 

ready for N s ?s 200 whereas CMc needs N s m 300 events 
to reach £ m m 0.71. Since the OMc method turns out to 
be more efficient that CMc in all the parameter ranges, 
henceforth we continue using only the OMc method for 
further estimations. 
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FIG. 3: Efficiency of OMc with a 40° angular cut for different 
values of Nb /N s . 

With neutron tagging, the value of Nb/N s decreases, 
and that increases the efficiency of the OMc method. In 
Fig. |3 we show the efficiency of OMc at different values 
of Nh/N s . At Nh/N s = 10, which corresponds to the 
tagging efficiency £ tag ~ 67%, the asymptotic efficiency 



of OMc already increases to 0.85, and at Nf,/N s = 5, 
corresponding to e tag ~ 83%, it reaches 0.90. Moreover, 
with decreasing N b /N s the asymptotic value is reached 
at lower and lower number of signal events. For higher 
tagging efficiencies, the optimal value of the angular cut 
increases. In fact, as noted before, in the limit of no back- 
ground the OM method without the cut is more optimal. 
However this limit is physically not reached due to the 
presence of oxygen events. 

The OMc method thus sacrifices less than 25%, and at 
higher tagging efficiency, even less than 15%, of the point- 
ing accuracy of the MLE method. On the other hand, 
it has the great advantage of being independent of the 
detailed neutrino spectra and cross sections. Therefore, 
this method can be extremely useful for a fast analysis of 
the SN signal. After all, an early warning would depend 
on a quick and simple data analysis while later one can 
certainly optimize by fitting detailed energy spectra to 
the observed signal. 

We stress that our preference for a parameter-free 
method over MLE in this analysis is strongly influenced 
by the current status of our knowledge regarding the pdf 
of the angular distribution. It may indeed be possible 
to use MLE and include all the systematic uncertainties, 
perhaps giving a better estimate for the pointing accu- 
racy. However, faced with a tradeoff between model in- 
dependence and higher efficiency, we give more weight to 
the former. If in future we understand the primary spec- 
tra much better than we do now, this preference may 
change. 



III. SUPERNOVA POINTING ACCURACY OF 
WATER CHERENKOV DETECTORS 

We now apply our method to a more realistic repre- 
sentation of the SN signal in a water Cherenkov detec- 
tor. We shall limit our analysis of the pointing accuracy 
to our best parameter-free method, i.e. the Orientation 
Matrix method with an angular cut (OMc). 

In order to determine the pointing accuracy numer- 
ically we simulate a large ensemble of SN signals in a 
water Cherenkov detector, assuming different efficiencies 
for neutron tagging. To this end we assume that the SN 
is at a distance D = 10 kpc and releases the neutron- 
star binding energy Ej, = 3 X 10 53 erg in the form of 
neutrinos. Details of the assumed neutrino spectra and 
fluxes are given in Appendix lAl The spread in the pre- 
dicted neutrino spectra has been taken care of by using 
two models, a model from the Garching group (model 
G) J23j and a model from the Livermore group (model 
L) [24| as described in the same appendix. We take into 
account the effects of neutrino flavor conversions by con- 
sidering the three mixing scenarios, (a) normal mass hier- 
archy and sin 2 13 > 10 -3 , (b) inverted mass hierarchy 
and sin 2 13 > 10~ 3 , and (c) any mass hierarchy and 
sin 2 9 13 < 10~ 3 . The six combinations of the models 
and neutrino mixing scenarios are represented by G-a, 
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G-b, G-c, L-a, L-b, L-c. We use sin 2 (26 ) = 0.9 for the 
solar neutrino mixing angle. 

As reaction channels we use elastic scattering on elec- 
trons ve~ — ► ver ', inverse beta decay D e p — ► ne + , and 
the charged-current reaction v e + 16 — ► X + e _ , while 
neglecting the other, subdominant reactions on oxygen. 
The cross sections for these reactions arc summarized in 
Appendix [BJ The oxygen reaction is included because 
it provides the dominant background for the directional 
electron scattering reaction in a detector configuration 
with neutron tagging where the inverse beta reaction can 
be rejected. 

For the detector we assume perfect efficiency above an 
"analysis threshold" of 7 MeV, and a vanishing efficiency 
below this energy. The actual detector threshold may be 
as low as 5 MeV. Though lowering the threshold increases 
the ratio of elastic scattering events and the inverse beta 
events, it also introduces a background from the neutral- 
current excitations of oxygen (see Appendix[BJ). In order 
to avoid additional uncertainties from the cross section of 
these oxygen reactions, we use the higher analysis thresh- 
old. We have checked that the net improvement by low- 
ering the threshold to 5 MeV is less than 10% in all cases. 

We assume a fiducial detector mass of 32 kiloton of wa- 
ter. Using the neutrino spectra and mixing parameters 
from the six cases mentioned above, we obtain 250-300 
electron scattering events, 7000-11500 inverse beta de- 
cays, and 150-800 oxygen events. The ranges correspond 
to the variation due to the six different combinations of 
neutrino mixing scenarios and models for the initial spec- 
tra. 

The procedure of event generation is described in Ap- 
pendix El The angular resolution function of the Super- 
Kamiokande detector does not follow a Gaussian distri- 
bution, rather it is close to a Landau distribution that we 
use for our simulation. In Fig.U the angular distribution 
of v e p — > ne + events (green) and elastic scattering events 
ve~ — ► ve~ (blue) of one of the simulated SNc arc shown 
in Hammer- Aitoff projection, which is an area preserving 
map from a sphere to a plane. 

The position of the SN is estimated with the OMc 
method. As explained in Sec.|Hj the optimal value of the 
angular cut depends on the neutron tagging efficiency as 
well as the neutrino spectra. We use a sharp cutoff with 
30° opening angle for the OMc, which may not be opti- 
mal, but is observed to be close to optimal in almost the 
whole parameter range. For low values of etag, the value 
of the cut should be lowered whereas for large values of 
Stag it should be increased by about 10°. The optimal 
cut depends also on the details of the detector properties 
and neutrino spectra. 

A histogram of the angular distances between the true 
and the estimated SN position found in 40000 simulated 
SNe for different neutron tagging efficiencies for the case 
G-a is shown in Fig. [5J The histogram fits well the dis- 
tribution 




FIG. 4: Angular distribution of v e p — 
and elastic scattering events ve~ — ► ve 
lated SN. 



ne + events (green) 
(blue) of one simu- 



wherc £ is the angle between the actual and the estimated 
SN direction, and S is a fit parameter. 




FIG. 5: Histogram of the angular distance £ of the esti- 
mated SN direction to the true one for 40000 simulated SNe 
with neutrino parameters corresponding to G-a and neutron 



tagging efficiencies eta 



0,0.8 and 1. The fits using the 



(15) 



distribution in Eq. 1151 are also shown. 

Defining the opening angle £ a for a given confidence 
level a as the value of I for which the SN direction esti- 
mated by a fraction a of all the experiments is contained 
within a cone of opening angle £, we show in Fig. |G| the 
opening angle for 95% C. L. for the six cases of neu- 
trino parameters. Clearly, the pointing accuracy depends 
weakly on the neutrino mixing scenario as well as the 
initial neutrino spectra. Some salient features of this de- 
pendence may be understood qualitatively as follows. 

The signal events are dominated by v e . Indeed, nearly 
half of the clastic scattering events are due to v e , whereas 
the remaining half are due to the other five neutrino 
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FIG. 6: The pointing accuracy ^95 as a function of the neu- 
tron tagging efficiency e ta g for six cases corresponding to three 
neutrino mixing scenarios and two models for the initial neu- 
trino spectra. 



species. The cross section of electron scattering events 
increases with energy. Therefore, the more energetic the 
v e arriving at the detector, the larger the number of sig- 
nal events and the better the pointing accuracy. Though 
the initial average v e energies are equal in the models G 
and L, the model L gives a much larger average energy 
for the initial v x spectrum. The v e -v x mixing then tends 
to give more energetic v e in the model L. As a result, 
for each mixing scenario, the model L predicts a better 
pointing accuracy than the model G. 

Within a model, the pointing accuracy is governed by 
the background-to-signal ratio. Since the cross section of 
the dominating background v e p reaction increases with 
energy, more v e -v x mixing tends to give more energetic 
v e and hence more background and less pointing accu- 
racy. The ratio of v e -v x mixing to v e -v x mixing within a 
model is the smallest for the mixing scenario (a) and the 
largest for the scenario (b). Therefore, within a model 
the scenario (a) always gives the best pointing accuracy 
and the scenario (b) always gives the worst. Note that in 
the limit etag = 1.0 when all the v e p background is elim- 
inated, the scenarios (b) and (c) give identical pointing 
accuracies since the final v e spectra in these two schemes 
are identical. 

If neutrino experiments or supernova simulations have 
already identified the actual scenario, we could use the 
confidence limits given by that particular scenario. In- 
deed, the information for identifying the mixing scenario 
may already be contained in the observed neutrino spec- 



TABLE I: Opening angle t a of the cone with a confidence 
level to contain the true SN direction for different tagging 
efficiencies with the "worst case" scenario. The bottom row 
gives the width 5 of the Gaussian distribution f(£). The * 
column gives the pointing accuracy in the limit when all the 
background, including the oxygen events, is weeded out. 
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0.5 


0.8 


0.9 


0.95 


1.0 
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= 0.68 
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3.0 
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1.7 





= 0.90 
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2.9 
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= 0.95 


7.8 
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3.2 


3.0 


2.9 
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= 0.99 


fO.O 


6.1 


4.5 


4.1 


3.9 


3.7 


3.6 
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3.0 


1.9 


1.4 


1.3 


1.2 


1.2 


1.1 



tra themselves which may be extracted by further data 
analysis SIS 00. However, in the absence of the 
immediate availability of this information, one has to se- 
lect the least efficient scenario, G-b, in order to obtain 
the most conservative limits. This is the scenario with 
the least v e -v x mixing, the largest v e -v x mixing and the 
lowest initial average energy for v x . 

Since the pointing accuracy is worse when the com- 
ponent of initial v x flux in the final v e flux is smaller, 
the "worst case" scenario is expected to be when the v e 
survival probability p (see Appendix^) takes its highest 
possible value, which is nearly 0.45 at 3tr [2f| , corre- 
sponding to sin 2 (26 ) ps 0.99. In Table [Q we give the 
values of the opening angle £ for this "worst case" sce- 
nario for various confidence levels a and tagging efficien- 
cies. We also give the values of the fit for 6 in Eq. (TSJ, 
from which the numbers for any confidence level can be 
read off. For e tag = 0, at 95% C.L. the pointing accuracy 
is 7.8°, which improves to 3.6° for e t ag = 80% and 3° for 
^ tag = 1 • This is nearly a factor of 3 improvement in the 
pointing angle, which corresponds to almost an order of 
magnitude improvement in the area of the sky in which 
the SN is located. 

The last column (marked by *) in the table shows the 
pointing accuracy in the limit of no background, i.e. the 
case where all the inverse beta decay as well as the oxygen 
events are weeded out. This gives the intrinsic limiting 
accuracy due to the angle of electron scattering, the an- 
gular resolution of the detector and the efficiency of our 
OMc algorithm. For 95% C.L. the table gives a pointing 
accuracy of 2.9°. This may be compared with the 2.4° 
that was estimated in Sec. Ill Al for the toy model when 
there was no background. The degradation in the point- 
ing accuracy may be attributed to the loss of information 
in the OMc method, the 10% smaller number of events in 
the SN simulation, and the difference between the actual 
angular distribution and the Gaussian as taken in the toy 
model. 

For a SN at 10 kpc, in the worst case scenario we get 
nearly 10600 events, out of which the electron scatter- 
ing signal is N s ss 270. We are then already at or just 
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below the asymptotic limit for the efficiency of the OMc 
method. For a larger detector like Hyper-Kamiokandc, 
the desired accuracy can be calculated simply by rescal- 
ing according to the number of signal events. For a 
detector with 25 times the fiducial volume of Super- 
Kamiokande, the pointing accuracy is then expected to 
be 2° without gadolinium and 0.6° with e t ag > 90%. If 
the total number of events is much smaller than 10000, 
for small £t a g the efficiency of OMc will be smaller and 
this factor needs to be taken into account for calculating 
the real pointing accuracy. For e tag > 0.8, the scaling 
with number of events should even work for a very small 
number of events, as can be seen from Fig. [3] 

Note that all the model dependences discussed here, 
though clearly observable, are only around 10% (see 
Fig. HJ). This indicates that the pointing accuracy es- 
timates are quite model independent, and hence robust. 



IV. HIGH-ENERGY NEUTRINO TELESCOPES 
A. High-energy supernova neutrinos 

Turning to the putative TeV neutrino burst associated 
with a SN explosion we note that the shock wave may 
well accelerate protons to energies up to 10 16 eV. This 
idea is supported by the fact that the cosmic rays be- 
low the knee, E < 10 16 eV, contain a total amount of 
energy comparable to that injected by all galactic SNe. 
The possibility of detecting high-energy neutrinos from a 
galactic SN has been discussed in the literature. In par- 
ticular, it was shown that during the first year after the 
explosion, the SN shock wave will produce a large flux of 
neutrinos with energies above 100 GeV, inducing more 
than 10 3 muons in a km 2 detector [2||. More recently it 
was suggested that the high-energy neutrino signal would 
arrive just 12 hours after the SN explosion and would last 
for about one hour, giving about 100 muon events with 
E > 1 TeV in a km 2 detector Q. 

Of course, the number of expected events strongly de- 
pends on unknown parameters, in particular on the total 
energy emitted in the form of pions £^ ot at a given time 
and the maximum energy of the emitted neutrinos E m&x . 
The neutrino spectrum depends on that of the acceler- 
ated protons. However, higher-energy neutrinos have a 
greater chance of being detected so that we are not in- 
terested in the exact spectral shape once it is a typical 
power-law for shock acceleration, dN v /dE cx E~ a with 
a < 2. If the spectrum is softer, a > 2, detecting the 
neutrinos is more difficult. Note that the spectrum of 
protons after shock acceleration can be dominated by 
high-energy particles in some cases |27^ . or even can be 
monochromatic if protons are accelerated in a potential 
gap. We will concentrate on the "standard" case of a 
E~ 2 neutrino spectrum, although our results for other 
cases will be similar. Following the calculation of Ref. 
and assuming -E max = 1 TeV we expect around 50 muon 
events in a km 2 detector during 1 hour at a time about 



12 hours after the SN explosion. For a larger maximum 
energy, E max = 1000 TeV, the number of muons increases 
to 200. 



B. Signal from below 

In order to discuss the expected signal in different neu- 
trino telescopes j^, we begin with the case where the 
SN happens in a part of the sky that a given neutrino tele- 
scope sees through the Earth. Future km 2 detectors like 
IceCube at the South Pole [lij, or northern projects like 
the Gigaton Water Detector at Baikal [22 and Nemo in 
the Mediterranean can detect the high energy neu- 
trinos. For each event the angular resolution is around 
one degree. In this case the pointing to the SN can be 
resolved with an accuracy of about l°/v / 50 « 8' (arc- 
minutes). However, this purely statistical error does not 
include possible systematic effects. Most important is 
the limited knowledge of the alignment of these detec- 
tors. Therefore, the pointing accuracy of a km 2 detector 
is probably larger and around a few tenths of a degree. 

Even the existing smaller detectors can see a significant 
signal. The northern sky is under control of AMANDA-II 
with an effective area of 0.1 km 2 and angular resolution of 
2° at TeV energies. AMANDA-II will then detect 5 (20) 
events for E max = 1 TeV (1000 TeV) and thus will be 
able to resolve the SN direction to better than 1°. After 
their completion, the northern projects ANTARES [34[ 
and NESTOR [3f| will be comparable to AMANDA-II. 



C. Signal from above 

If the high-energy SN neutrinos arrive "from above," 
they are masked by the large background of at- 
mospheric muons. For IceCube this background is 
about 5.2 x 10 10 year -1 from the upper hemisphere at 
the "trigger" level |3(|. This corresponds to nearly 
300 hour -1 degree -2 . If we note that the angular res- 
olution is about 1°, the expected signal of 100 events 
will be much larger than the background fluctuations in 
one pixel of the sky. Morever, the expected SN neutri- 
nos will have multi-TeV energies so that energy cuts will 
reduce the background. Also, the significance of the sig- 
nal will be enhanced by the angular prior defined by the 
low-energy signal in Supcr-Kamiokande. 

For AMANDA-II size detectors both background and 
signal are about 10 times smaller. Moreover, the angu- 
lar resolution is only about 2°. Therefore, the expected 
signal in one pixel of the sky will be comparable to the 
background fluctuations. Energy cuts may allow one to 
detect the SN signal from the "bad" side of the sky. 
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V. SUMMARY AND CONCLUSIONS 

The MeV neutrinos from the cooling phase of a SN 
will arrive at the Earth several hours before the optical 
explosion. These neutrinos will not only give an early 
warning of the advent of a SN explosion, but they can 
also be used to determine the location of the SN in the 
sky, so that the optical telescopes may concentrate on a 
small area for the observation. 

In a water Cherenkov detector like Supcr-Kamiokandc, 
the v-e scattering events arc forward peaked and thus can 
be used for the pointing. The main background comes 
from the inverse beta decay reactions v e p — ► ne + , which 
is nearly isotropic and has a strength more than 30 times 
that of the electron scattering "signal" reaction. The 
reactions of the neutrinos with oxygen also contribute to 
the nearly isotropic background. 

The authors of Ref. [Tl| have estimated the pointing 
accuracy using the Rao-Cramer bound. This may over- 
estimate the accuracy by ~ 10% since even the most 
efficient method, the maximum likelihood estimate, can 
reach the minimum variance bound only for a large num- 
ber of signal events. More importantly, the maximum 
likelihood method needs as an input the exact form of 
the fit function, which is not available due to our cur- 
rently poor knowledge of the neutrino spectra. Taking 
into account the errors due to a wrong choice of the fit 
function is difficult. Therefore we choose to calculate the 
pointing accuracy using a concrete and simple estima- 
tion method that is independent of the form of the fit 
function. 

We explore some parameter-free methods that only use 
the data and exploit the symmetries inherent in the phys- 
ical situation, and therefore give a model independent es- 
timation of the pointing accuracy while sacrificing some 
information from the data. We perform a statistical anal- 
ysis of these methods using a toy model, which has a 
Gaussian signal on top of an isotropic background. We 
find that a method that uses the "orientation matrix" 
with an appropriate angular cut (OMc) is an efficient 
method that uses more than 75% of the information con- 
tained in the data if N s > 200. We argue that this loss 
of information is well worth the gain of model indepen- 
dence, and continue to use this method for the simulation 
of the actual angular distribution at a water Cherenkov 
detector. It turns out that this method is much more 
efficient than the one used in Ref. (l2T| . 

One may add gadolinium to Super-Kamiokande in or- 
der to tag neutrons and therefore reduce the background 
due to inverse beta decay events. We quantify the in- 
crease in the pointing accuracy as a function of the neu- 
tron tagging efficiency £t ag . It is found that the accuracy 
increases by more than a factor of two with £ tag = 0.8 and 
by nearly a factor of 3 for e tag = 0.95. For e tag > 0.95, 
the oxygen events act as the major background and that 
saturates the advantage of increasing the tagging effi- 
ciency beyond this value. 

The efficiency of the OMc method improves with a 



smaller background-to-signal ratio, which makes this 
method even more useful at high tagging efficiencies. It 
is also observed that at higher e tag this method attains 
its maximum efficiency level for a much smaller number 
of events. The optimal value of the cut increases with 
increasing £t a g, though this dependence is weak and we 
perform our estimations with a fixed value of the cut, 
which is near optimal in the whole parameter range. Our 
estimations are therefore slightly conservative. 

With a simulation of a water Cherenkov detector like 
Super-Kamiokande, we determine the pointing accuracy 
obtained from a SN at 10 kpc. The accuracy has a weak 
dependence on the neutrino mixing scenarios and the ini- 
tial neutrino spectra. We find that the worst case sce- 
nario is when the detected 9 e spectrum has the largest 
admixture of the initial v x spectrum, and the detected 
v e spectrum has the lowest admixture of the initial v x 
spectrum. This worst case turns out to be the one with 
the inverted neutrino mass hierarchy, sin 2 13 > 10~ 3 , 
and the largest possible solar mixing angle. The OMc 
method gives the pointing accuracy of 7.8° at 95% C.L. 
without neutron tagging. The accuracy improves to 3.6° 
at 80% tagging efficiency and to 3.2° at 95% tagging ef- 
ficiency. Beyond this, the pointing accuracy saturates 
due to the presence of the oxygen events and the limited 
angular resolution of the detector. For a larger detec- 
tor, the expected accuracy may be scaled according to 
the number of events. At a megaton detector like Hyper- 
Kamiokande, this gives an accuracy of 0.6° for e ta g > 0.9. 

The SN shock wave may produce a TeV neutrino burst 
that arrives at the Earth within a day of the initial MeV 
neutrino signal. This can give about 100 events with 
E > 1 TeV at a km 2 detector like IceCube. Since the 
angular resolution of this detector is as good as 1°, the SN 
may be located to an accuracy of a few tenths of a degree. 
The limiting factor here is the alignment error of these 
detectors. The time correlation and the directionality 
of the events allows IceCube to detect them even "from 
above" against the background of atmospheric muons. 
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TABLE II: The parameters in the neutrino spectra models 
from the Garching group and the Livermore group. 



Model 


<**>(".)) 


(E {Ve)) 




$o(fe) 


$o(^e) 
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12 


15 


18 


0.8 


0.8 


L 


12 


15 


21 


2.0 


1.6 



TABLE III: The possible combinations of survival probabili- 
ties p and p. 

Case Hierarchy sin 2 O13 p p 

(a) Normal > 10" 3 cos 2 9 Q 

(b) Inverted > 10~ 3 sin 2 9 

(c) Any < 10" 3 sin 2 9© cos 2 6 Q 



APPENDIX A: NEUTRINO FLUXES 

For the time-integrated neutrino fluxes we assume dis- 
tributions of the form (33 



F° 



$0 (1 

E r(i + a) 



E 



exp 



-(a + 1 



E 



(Al) 

where F° denotes the flux of a neutrino species emitted 
by the SN scaled appropriately to the distance travelled 
from the SN to Earth. Here Eq is the average energy and 
a a parameter that relates to the width of the spectrum 
and typically takes on values 2.5-5, depending on the 
flavor and the phase of neutrino emission. The values of 
the total flux $0 an d the spectral parameters a and Eq 
are generally different for v e , v e and v x , where v x stands 
for any of v^ r or v^t- 

We consider two models for the initial neutrino fluxes. 
The first one is the recent calculation from the Garch- 
ing group j23, which we refer to as the model G. It 
includes all relevant neutrino interaction rates, includ- 
ing nucleoli bremsstahlung, neutrino pair processes, weak 
magnetism, nuclcon recoils and nuclear correlation ef- 
fects. The second one is the Livermore simulation [24|, 
referred to as model L, that represents traditional pre- 
dictions for flavor-dependent SN neutrino spectra that 
have been used in many previous analyses. The param- 
eters of these models are shown in Tabic |n] We take 
a(v e ) = a(P e ) = a(Vx) = 3.0 for both models. The value 
of a is not expected to have any significant influence on 
the pointing accuracy. 

When neutrino mixing is taken into account, the fluxes 
arriving at a detector are 



= pFl +(l-p)F° 



(A2) 
(A3) 

(1 - + (1 - p)Fl + (2 + p + p)Fl . (A4) 



Since the four neutrino species v x cannot be distinguished 
at the detectors, we only give the sum of their fluxes, 
AFv^ . Here p and p are the survival probabilities of v e 
and v e respectively. 

Depending on the mass hierarchy and the value of the 
mixing angle 613, the survival probabilities p and p be- 
long to one of the three mixing scenarios shown in Ta- 
ble IIIII We neglect the Earth matter effects and the 
details of the "transition" region around sin 2 813 ~ 10 -3 
|3EI- We also neglect terms of order (013) 2 in p and p. 



TABLE IV: Coefficients used in Eq. iBll for the elastic scat- 
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Then p and p depend only on the solar mixing angle as 
given in the table, and are independent of the values of 
solar and atmospheric mass squared differences as well 
as the atmospheric mixing angle. 



APPENDIX B: 



NEUTRINO REACTIONS IN 
WATER 



1. Elastic Scattering on Electrons 

The differential cross-section of the reaction v + e" 
v + e~ with v — {v e , v e , ^ jT , ^. r } is given by 



da G 2 F m e E v 



dy 



2tt 



A + B(l-y) 



C E V 



(Bl) 



where y = E e / E v is the energy fraction transfered to the 
electron, Gf the Fermi constant, and m e the electron 
mass. The coefficients A, B and C differ for the four 
different reaction channels and are given in Table 
The vector and axial-vector coupling constants have the 



2 sin 0\v and Ca 



with 



usual values Cv = - 
sin 2 Q w « 0.231. 

The scattering angle $ between the incoming neutrino 
and final-state electron is implied by 



2 {m e /E v ) cos 2 d 
(l + m e /E v ) 2 - cos 2 ?? 



(B2) 



2. Inverse Beta Decay 

For v e p — > ne + we use the differential and total cross- 
sections Eqs. (5-7) of Ref. [33, including the leading 
QED radiative corrections. 
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3. Oxygen as a Target 

Another important reaction is the charged-current v e 
absorption on oxygen |39l l4p| | . The dominant channels 
are 



V e ~ 


- 16 - 


15 0+p + e~, 


(B3) 


Ve ~ 


~ 16 - 


■+ 15 0* +p + 7 + e~, 


(B4) 


V & ~ 


- 16 - 


■* u N*+ P + p + e-. 


(B5) 



While these reactions cause far fewer events in a wa- 
ter Cherenkov detector than inverse beta decay, they do 
not have final-state neutrons and thus cannot be tagged. 
Therefore, in a detector configuration with efficient neu- 
tron tagging, these reactions provide the dominant back- 
ground to the directional electron scattering reactions. 

The neutrino energy threshold in these reactions is ap- 
proximately 15 MeV. The total cross section, summed 
over all channels, has been tabulated for the range 
15 < E v < 100 MeV Hjj. Directly above thresh- 
old the cross section is very small. We find that for 
25 < E v < 100 MeV the tabulated cross sections are 
nicely represented by the analytic fit 



16 



4.7 x 10" 40 cm 



MeV 



15 



1/4 



(B6) 



For an accurate determination of the detector response 
one needs the differential distribution of final-state e~ 
energies and angular directions for a given incident E v . 
Ref. [3^ provides extensive plots of such distributions 
after folding them with thermal E v distributions. This 
information is too indirect for our purposes. Therefore, 
we limit our investigation to a schematic implementation 
of this process where we assume that in every reaction the 
final-state energy is E e = E v — 15 MeV. For the angular 
distribution we assume 



da 



dcosfi 



1 



1 + (E e /25 MeV) 4 
3 + (Se/25 MeV) 4 



COS l? . 



(B7) 



where $ is the angle between incident v e and final-state 
e~ . This means that for small energies the angular distri- 
bution is proportional to 1 — | cos $ while for large ener- 
gies it is 1 — cos?9, i.e. it becomes more backward peaked 
for high energies. Our schematic approach roughly mim- 
ics the behavior shown in Ref. Since the oxygen 
cross section is very energy dependent, the contribution 
of this reaction to the pointing accuracy depends sensi- 
tively on the neutrino energy spectrum and is thus very 
uncertain anyway. 

Another potentially important class of charged-current 
reactions is HM l40l 



16 



O^X 



(B8) 



However, the contribution to the detector signal is some- 
what smaller than caused by the above v e reactions. 



Moreover, the v e reactions typically involve final-state 
neutrons and thus are rejected by neutron tagging. One 
exception is 



'O 



5 N- 



(B9) 



but its contribution is small. Therefore, we neglect this 
entire class of reactions in our study. 

Another class of reactions is the neutral-current exci- 
tation of oxygen |4lj 



3 -> v + X + 7 



(BIO) 



Most of these reactions cannot be rejected by neutron 
tagging. However, the total cross section for neutral- 
current scattering, including the channels without final- 
state 7, is smaller than for the charged-current v e reac- 
tion |4J|. Moreover, the 7 energies are below 10 MeV, 
and most of them even below our analysis threshold of 
7 MeV. Therefore, we also neglect this class of reactions. 



APPENDIX C: EVENT GENERATION 

For the generation of each of the events the following 
steps are performed: 

1. The energy E v of the reacting neutrino is chosen 
according to F v (E v )a{E v ). The energy E e and the 
scattering angle of the outgoing electron/positron 
is chosen according to the differential cross-section 
of the particular reaction. 

2. The measured energy E^t of the scattered elec- 
tron/positron is determined by adding Gaussian 
noise with variance a(E e ) = \/E £ Eq, where Eq = 
0.22 MeV. If E dct < E t h = 7 MeV, the event is not 
used in the data analysis. 

3. The measured position (0, i?) of the event is simu- 
lated according the angular resolution function of 
S up er- K amiokande . 

The angular resolution lZ(£)d£ is defined as the prob- 
ability that inside two cones with opening angles £ and 
£ + dl around the true direction the reconstructed direc- 
tion is contained. Reference |2lj gives numerical values 
for the opening angle £ of a cone around the true direc- 
tion which contains 68% of the reconstructed directions 
as well as the values of TZ(£) as a function of £ for var- 
ious energies. An inspection by eye shows that 1Z(£) is 
characterized by a large tail and cannot be well fitted by 
a Gaussian distribution. Inspired by the Landau distri- 
bution for energy losses, we have found that TZ(£) is well 
described by 



Tl{l) =Cexp - 



sin(^) 



(CI) 



12 

1 



12 
10 

8 
6 
4 
2 




where C is a normalization constant, and 



M E=16.1MeV 
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FIG. 7: Our fit of the angular resolution 1Z(£) as a function 
of I for positron energies E e — 6.8 and 16.1 MeV together 
with measured points from Ref. l2l| . 
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For illustration, wc show in Fig. [7] the angular resolution 
as a function of I for positron energies E e = 6.8 
and 16.1 MeV as implemented in our simulation, together 
with the measurements extracted from El. 
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